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Abstract 

_  J 

The  crushing  analysis  of  rotationally  symmetric  plastic 
shells  undergoing  very  large  deflections  is  presented.  A 
general  methodology  is  developed  and  simple  closed  form  solu¬ 
tions  which  can  be  useful  for  practical  applications  are  derived 
for  the  case  of  a  conical  shell  and  a  spherical  shell  under 
point  load,  a  spherical  shell  crushed  between  rigid  plates  and 
under  boss  loading,  and  a  spherical  cap  under  external  uniform 
pressure.  The  effect  of  the  end  conditions  and  the  limitations 
of  this  approach  are  discussed  in  detail. 
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1 .  Introduction 

Thin  metal  shells  of  revolution  subjected  to  axisyiranetric 
compressive  forces  are  capable  of  carrying  substantial  loads 
when  deflections  exceed  by  one  or  two  orders  of  magnitude  the 
thickness  of  the  shell,  and  become  comparable  to  the  largest 
linear  dimension  of  the  structure.  Thus,  the  deformation 
process  is  far  beyond  what  is  commonly  understood  as  post- 
buckling  behavior  of  elastic  or  plastic  shells.  A  distinctive 
feature  of  such  deformation  mechanism  is  that  the  strain  energy 
function  (in  the  case  of  elastic  shells) ,  or  dissipation  function 
(in  the  case  of  plastic  shells) ,  is  concentrated  in  narrow  zones 
(fold  lines  or  hinge  lines)  while  the  remainder  of  the  struc¬ 
ture  is  undergoing  a  rigid  body  motion.  Furthermore,  with  some  de¬ 
cree  of  symmetry,  the  fold  lines  are  forced  to  move  as  the  deforma¬ 
tion  process  goes  on.  Depending  on  the  geometrical  parameters 
of  the  structure  and  imperfections,  either  a  symmetric  or  un- 
symmetric  deformation  pattern  may  develop,  even  though  the  loading 
and  the  structure  itself  are  axially  symmetric. 

A  known  example  of  non-symmetric  but  regular  deflection  field  is 
the  Yoshinura  buckle  pattern  in  conpressed  thin  elastic  shells.  In 
the  range  of  very  large  deflections,  the  geometry  of  such  a  field 
can  be  conveniently  described  using  the  concept  of  isometric 
transformation  of  surfaces.  The  foundations  of  the  respective 
mathematical  theory  were  laid  by  Pogorielov,  [1],  who  also 
presented  brilliant  applications  of  his  theory  to  the  analysis 


of  post-critical  behavior  of  thin  elastic  shells,  [2],  His 
work  was  later  followed  by  Lukasiewicz  and  Szyszkowski  [3], 
while  Foster  [4]  presented  a  parallel  but  independent  study 
of  the  compression  of  thin  tubes.  Pugsley  and  Macaulay 
appear  to  be  the  first  to  study  the  crumpling  process  of 
thin  cylindrical  shells  with  diamond  shape  lobes  in  the  plastic 
range  [5] .  Although  no  provision  was  made  to  accomodate  in 
the  theory  the  travelling  hinges,  the  corresponding  calcula¬ 
tions  for  the  mean  crushing  force,  especially  the  more  recent 
ones  [6],  agree  well  with  experimental  results.  Further  contri¬ 
bution  to  the  understanding  of  the  process  of  progressive 
crumpling  of  tubes  in  the  multi-lobe  modes  was  made  in  [7] 
and  [8],  but  no  satisfactory  solution  of  this  interesting  problem 
was  offered  to  date. 

As  tubes  get  thicker  (R/h  <  50)  ,  a  well-known  transition 
takes  place  from  the  diamond  pattern  to  the  crinkling  (concertina) 
mode  of  deformation,  discussed  by  Alexander  [7],  Assuming  a 
simple  deformation  mechanism  with  circumferential  extension 
and  equating  the  rate  of  internal  energy  dissipation  to  the 
rate  of  work  of  external  forces,  he  was  able  to  calculate  an 
average  crushing  force.  Minimizing  next  the  force  level  he 
evaluated  the  length  of  the  local  buckling  wave,  and  derived  a 
simple  expression  for  the  optimal  value  of  the  crushing  force. 

The  above  procedure  has  proved  to  be  most  effective  in  the 
approximate  analysis  of  the  crumpling  process  of  various  thin- 
walled  structures  [8,10]. 


The  tube  inversion  (outside-in  or  inside-out)  is  another 
example  of  an  axially-symmetric  deformation  of  tubes.  Two 
solutions  for  the  magnitude  of  axial  loading  necessary  to  main¬ 
tain  the  plastic  flow  are  available.  One,  due  to  Al  Hassani 
et  al.  [10]  ,  is  good  as  long  as  the  assumption  of  a  constant 
thickness  holds,  and  another  one,  due  to  Abraraowicz  [11],  takes 
into  account  changes  of  thickness  and  circumferential  curvature. 
The  problem  of  hemispherical  shell  or  spherical  cap  loaded  by 
central  point-load,  rigid  boss,  or  crushed  between  plates,  received 
much  less  attention  in  the  literature  than  the  related  buckling 
problem  under  external  pressure  loading.  Updike  and  Kalnins 
performed  a  thorough  analysis  of  postbuckling  behaviour  of 
elastic  spherical  shells  compressed  between  rigid  plates,  [12, 

13]  .  Using  F.eissner  equations  for  shallow  rotationally  symmetric 
shells,  they  determined  the  force-deflection  relationship  in 
subsequent  stages  of  deformation,  and  the  point  of  bifurcation 
into  a  symmetric  inward  dimple,  and  later  into  an  unsymmetric 
shape.  In  the  subsequent  study  [14]  Updike  presented  an  approxi¬ 
mate  solution  for  the  rigid-plastic  shallow  shell  undergoing 
symmetric  deformations.  An  extensive  theoretical  and  experi¬ 
mental  investigation  into  this  problem  was  made  by  Kitching 
et  al.  in  [15],  where  the  transition  into  an  unsymmetric  plastic 
mode  was  calculated  numerically.  Experiments  in  which  polygon 
rather  than  circular  circumference  dimple  were  observed,  were 
reportedly  made  earlier  by  Pian  [9]. 


Leckie  and  Penny  [16]  performed  a  series  of  tests  on 
carefully  manufactured  hemispherical  shells  loaded  centrally 
by  a  rigid  boss.  These  experiments  were  followed  by  a  theoretical 
study  due  to  Morris  and  Calladine  [17]  ,  in  which  the  deformed 
shape  of  the  shell  and  the  force-displacement  characteristics 
were  obtained  by  a  sequence  of  upper  bound  calculations  ,  in 
conjunction  with  an  ingeneous  idea  of  treating  the  shell  as  a 
three-dimensional  body,  [18] . 

The  work  of  Morris  and  Calladine  consititutes  a  milestone 
in  the  understanding  of  the  crushing  behaviour  of  shells  of 
revolution.  It  shows  through  relatively  simple  numerical  cal¬ 
culations  that  plastic  deformations  are  indeed  confined  to  a 
relatively  narrow  ring  or  section  of  a  toroidal  surface,  and 
that  the  ring  is  moving  outward  as  the  deformation  process 
goes  on,  leaving  a  rigid  region  behind.  A  snap-through  effect 
was  observed  whose  extent  depends  on  the  boss  size  p. 

The  plastic  behaviour  of  shallow  spherical  shells  under 
uniform  pressure  loading  is  similar  to  that  of  a  shell  under 
rigid  boss  loading,  as  far  as  the  occurrence  of  the  plastic 
snap-through  phenomenon  is  concerned.  This  problem  was  studied 
by  Puszek  [19]  and  Jones  and  Ich  [20],  using  an  approximate  yield 
condition  for  rotationally  symmetric  shells.  Recently,  the 
same  problem  was  reconsidered  by  Kondo  and  Pian  [21]  by  applyine 
the  generalized  yield  line  method,  v/hich  takes  into  account 
the  changes  in  the  geometry  of  the  structure . 


In  the  present  paper  we  shall  extend  the  approach  developed 
in  [7]  and  [17]  to  general  shells  of  revolution,  with  the 
purpose  of  deriving  some  simple  results  useful  for  practical 
purposes.  A  detailed  discussion  will  be  presented  on  the 
restrictions  imposed  on  the  solution  by  the  conditions  of 
kinematic  continuity  on  the  moving  hinge  circles,  and  inextensi 
bility  of  the  shell  in  the  meridional  direction..  Introducing 
simple  displacement  and  velocity  fields,  closed-form  solutions 
on  the  crushing  force  will  be  derived  for  a  variety  of  struc¬ 
tures,  including  conical  and  hemispherical  shells  compressed 
between  rigid  plates,  and  a  spherical  shell  under  a  point-load 
or  central  boss  loading.  The  influence  of  the  boss  size, 

and  clamped  edge  conditions  will  also  be  discussed  at 
length . 

The  analysis  will  be  restricted  to  axisymmetric  deforma¬ 
tion  modes  only.  A  simplified  approach  to  study  the  lateral 
crushing  of  cylindrical  shells  using  a  polygon  shaped  fan  of 
flat  elements  was  developed  by  Morris  and  Calladine  [ 22] . 

This  method  is  particularly  suited  for  treating  nonsymmetric 
but  regular  deformation  modes  in  hemispherical  shells  such 
as  the  ones  observed  by  Pian  [  8  ] .  The  treatment  of  non¬ 
symmetric  problems  with  smooth  deformation  fields  will  be 
the  subject  of  a  future  publication  [23]. 


2 .  Analysis  of  Discontinuities 

The  condition  of  kinematic  continuity  on  the  moving  hinges 
or  hinge  circles  in  rigid-plastic  beams  and  flat  circular  plates 
were  formulated  by  Hopkins  [24]  and  Hopkins  and  Prager  [25]. 

Before  writing  the  continuity  equation  for  a  general  shell 
consider  first  the  one-dimensional  case,  since  the  physical 
arguments  leading  to  these  equations  are  the  same. 

Let  T  denote  a  straight  line  yield  hinge  whose  position 
on  the  flat  surface  element  is  described  by  a  scalar  equation 
x  =  v(t).  The  downward  deflection  and  velocity  of  the  material 
point  lying  on  the  shell  middle  surface  are  denoted  respectively 
by  f(x,t)  and  f(x,t).  We  assume  that  deflections  and  velocities 
are  continous  at  any  time  and  at  any  point  of  the  shell. 

Since  the  deflection  f  is  continuous  across  F ,  the  time  deriva¬ 
tive  of  f(v(t),t)  along  T  is  also  continuous  across  F.  Thus 

[f]  +  0  (t)  [f '  ]  =  o  (i) 

where  the  brackets  denote  discontinuties ,  if  any,  of  the  enclosed 
quantities  across  r ,  i.e.  [n]  -  n+  -  n_.  Now,  f  is  continuous  across 
T  and  Eq.  (1)  therefore  shows  that  the  slope  f'  can  be  discontinuous 
only  across  stationary  yield  hinges,  for  which  v(t)  =0.  If  the 
yield  hinge  T  is  not  stationary,  then  f’  must  be  continuous  across 
T.  By  repeating  similar  arguments  one  can  show  that 

[f'l  +  v(t)  [ f ” ]  =  0  (2) 

provided  v  (t)  ^  0. 


Consider  now  a  hinge  line  moving  with  a  constant  velocity 


v(t)  =  V  down  the  undeformed  flat  surface  element 

leaving  a  plastically  deformed  region,  Fig.  1.  The  rate  of 

rotation  and  curvatures  on  both  sides  of  the  moving  hinge  are 

respectively  f+  =  0,  £'_  =  0,  f+  =  k,  f"  =  0.  Now,  Eq.  (2) 
reduces  to 

0  +  VK  =  0  (3) 

tve  can  see  that  the  moving  hinge  line  imposes  a  constant 
curvature  <  =  4  on  the  initially  flat  surface  element  if  the 
slope  of  the  velocity  field  6  is  related  to  the  hinge  velocity 
by  (3T\ 

Return  now  to  the  general  case,  following  the  derivation 
presented  in  [11]  and  consider  a  curved  discontinuity  line  F 
moving  down  the  middle  surface  of  the  shell  x*  The  line  r 
is  dividing  the  surface  x  into  two  separable  parts 
X  =  (  +  )  XUrU(  ^x*  Then  ^X  denotes  the  deformed  part,  subjected 
to  plastic  deformation.  The  vector  field  f  is  defined  over  the 
entire  surface  x-  Now,  the  previous  analysis  still  applies 
except  that  the  first  and  second  derivation  of  f  in  (1)  and 
(2)  should  be  replaced  by  the  evalutation  of  the  first  and 
second  gradient  of  f  on  the  surface  x.  This  is  equivalent 
to  the  replacement  of  partial  derivatives  in  (1)  and  (2)  hy 
covariant  derivatives,  so  these  eauations  now  read 


Tf1!  +  v  f 

fi|  1 

L  J  n  [ 

.  aJ 

(4) 


0 


i  =  U,2,3} 


-8- 


where  n  is  a  surface  vector  normal  to  the  discontinuity  line 
r,  while  Vn  denotes  the  velocity  of  the  hinge  line  in  the  n 
direction. 

Considerable  simplification  can  he  obtained  in  the  case 
of  rotationally  symmetric  shells.  The  hinge  line  becomes  now 
a  circle  and  each  of  the  continuity  condition  (4)  and  (5) 
yields  only  one  scalar  equation  in  the  meridional  plane 

[f]  +  v  ff,  ]  =  0  (C) 

n  a 


[f,a] 


v  [f. 

n  act 


where  [f,  ]  is  a  jump  in  the  rate  of  rotation  and  [f,  1  is 

a  -it-  aa 

the  corresponding  jump  in  the  principle  (meridional)  curvature. 

Take  now  an  intermediate  stage  of  the  crushing  process 
’•-'here  plastic  deformations  are  confined  to  a  relatively  narrow 
zone  contained  between  an  inner  and  outer  hinge  circle  B  and  C,  Fig. 
The  outer  hinge  circle  B  moves  down  the  undeformed  shell  with  a 
velocity  V  imposing  a  larger  curvature,  i.e.  changing  the 
local  radius  of  curvature  from  R_  to  rF.  The  material  in  front 

d  b 

of  the  hinge  is  rigid  and  undergoes  a  constant  and  finite  rate  or. 
rotation  0.  Thus,  the  continuity  condition  reduces  in  this 


case  to 

~“b  +  v44 '  h) = 


The  inner  hinge  C  is  travelling  with  the  velocity  and  remover 

the  curvature,  or  strickly  speaking,  it  changes  the  small  radius 

of  curvature  r  back  to  the  large  one  K  .  be  can  then  write 
C  C 


a  similar  equation  to  (8)  expressing  the  continuity  constraint 

u  +  V  I  —  -  =  0  ( 

c  c(r  c  bc; 

For  the  material  behind  the  circle  C  to  be  rigid,  the  shell 
should  translate  rather  than  rotate.  Hence,  there  must  be  at 
C  a  counterrotation  ^  of  the  same  magnitude  as  at  B.  The 
condition  ojg=oc  "‘implies  the  following  relation  between  the  six 
kinematic  and  geometric  parameters  of  the  problem 


Equation  (10)  can  be  further  simplified.  If  for  example  the 

deflected  shape  of  the  shell  is  assumed  to  be  isometric  to  the 

original  one,  it  must  be  its  mirror  reflection  i.e.  Rg  =  Rc 

see  [1].  Furthermore  the  shape  of  the  shell  middle 

surface  between  the  hinge  circle  B  and  C  may  be  approximated  by 

a  section  of  a  circular  toroidal  surface  so  that  r_  =  r  =  r, 

B  c 

as  suggested  in  [14,11]. 

Equation  (6)  carries  the  information  that  the  shape  of  shell 
at  B  and  C  should  be  continuous  at  all  times,  except  at  the 
stationary  hinge,  i.e.  at  the  rigid  boss  or  clamped  edge. 

In  the  next  section  we  shall  see  whether  conditions  (6) 
and  (10)  can  be  made  compatible  with  other  simplifying 
assumptions  regarding  the  kinematics  of  the  crushing  process. 

3.  Mechanisms  of  Plastic  Deformation 


The  foregoing  analysis  has  shown  that  there  are  two  con¬ 
centrated  plastic  hinges  at  B  and  C  and  that  the  cross  section 


-10- 


be  tween  them  rotates  as  a  rigid  shape  with  an  angular  velocity  w , 
as  the  displacement  of  the  shell  increases.  Our  aim  is  to  evaluate 
the  internal  energy  dissipated  in  such  a.  mode  of  deformation. 
Equating  this  energy  to  the  rate  of  work  of  external  forces 
an  instantaneous  crushing  force  can  be  found  necessary  to 
maintain  the  plastic  flow. 

The  dissipation  is  due  to  a  discontinuous  velocity  field 
at  plastic  hinges,  and  to  a  continuous  deformation  field  in  the 
shell  section  between  B  and  C.  Consider  first  the  latter.. 

In  rotationally  symmetric  shells  there  are  four  components 
of  the  generalized  strain  rate  and  stress  fields.  The  meridional 

curvature  rate  <  is  infinite  at  hinge  circles  and  vanishes 

•D 

between  them.  This  follows  from  the  assumption  that  the  rate 
of  rotation  is  a  step  function. 

It  should  be  noted  that  the  meridional  curvature  itself 
<5,  governed  by  (8)  is  kept  constant  at  any  material  point 
passing  through  the  plastically  deforming  annulus,  but  does  net 


have  to  be  constant  along  spatial  coordinate  between  B  and  C. 


The  circumferential  curvature  rate  <e  is  continuously 

changing  between  B  and  C  from  positive  —  through  zero  to  negative 
i_  r'B 

Rr/ and  consequently  some  energy  is  dissipated  in  this  deformation 


mechanism.  An  exact  formula  for  kp  was  derived  in  [11]  using 
a  rigorous  Eulerian  description.  In  the  same  paper  it  was 


shown  that  the  associated  energy  dissipation  is  a  small  fraction 


of  the  total  energy. 
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For  the  sake  of  simplicity  we  shall  disregard  <  in  the 

9 

equation  for  the  energy  dissipation-  A  very  good  estimate  of 
the  contribution  of  ic  to  that  energy  can  be  given  in  a 
straightforward  manner. 

Suppose  the  circumferential  curvature  is  reversed  between 
B  and  C  from  4  to  .  The  rate  of  work  done  is  twice 
the  energy  required  to  remove  the  curvature  from  4  to  zero . 

The  curvature  is  removed  by  rotating  two  ends  of  a  cut  circle 
through  the  angle  2 rr ,  (cf.  [26]).  Thus,  the  rate  of  energy 
in  the  considered  mechanism  of  deformation  can  be  approximated 
by  : 

E(k0)  =  2  (2tt)  MqV  =  4ttMq  |  w  (11) 

where  it  was  assumed  that  there  are  no  interaction  between  the  cir¬ 
cumferential  bending  moment  and  the  remaining  non-vanishing 
components  of  the  generalized  stress  field.  The  velocity  V  in 

(11)  is  an  average  downward  velocity  of  the  deforming  zone, 
b  • 

V  =  j  to;  see  Fig .  3 . 

At  certain  stages  of  the  crushing  process  considerable 
meridional  membrane  forces  may  be  developed.  In  the  present 
paper  the  meridional  extension  rate  is  assumed  to  vanish  in 
the  zone  of  continuous  deformations.  Any  extension  or 
shortening  of  the  arc  length  in  the  meridional  direction  is 
accommodated  by  the  hinge  circle  B  and  C.  Thus,  we  do  not 
require  inextensibility  of  the  material  in  the  meridional 
direction  but  instead  extensional  plastic  hinges  are  intro¬ 
duced.  This  assumption  is  in  the  spirit  of  Calladine' s:  approxi¬ 
mation  in  the  problem  of  moderately  large  deflection  of  clamped 
circular  plates  [18] . 
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The  only  significant  component  of  the  generalized  strain 
rate  is  the  circumferential  compression  A  .  In  order  to  calcu¬ 
late  the  associate  dissipation  energy  v;e  shall  strictly  follow 
the  approach  developed  in  [17],  and  treat  the  shell  section 
between  B  and  C  as  a  three-dimensional  tody.  The  rate  of  internal 
energy  dissipation  is  equal  to 


E  = 


/Vedv 

V 


a 

o 


x  e  9  dS 


where  V  and  S  denote  respectively  the  total  volume  and  area 
of  the  cross  section  between  B  and  C,  and  x  denotes-  a  current 
distance  of  the  given  point  F  from  the  axis  of  symmetry  of  the 
shell,  see  Fig.  3.  Note,  that  the  yield  stress  is  positive  in 
tension  zones  and  negative  in  compression  zones.  The  circumfer- 
ential  strain  rate  is  given  by  Eg  =  where  y  is  the  current 
distance  from  the  axis  of  instantaneous  rotation  I-I .  The 
position  of  this  axis  with  respect  to  the  mell  cross  section 
is  yet  to  be  determined.  Substituting  ea  in  (12)  yields 


E.  . 
int 


=  2710 


'o-/l 


y  |cS 


(13) 


One  has  to  distinguish  now  two  special  caeec.  If  the  rise  of 
the  toroidal  surface  over  the  points  L  anc.  C  is  comparable 
to  the  shell  thickness,  the  assumption  about  non-uniform 
distribution  of  ta  across  the  thickness  is  justified,  and  (13) 
holds.  If  the  converse  is  true,  one  can  replace  dS  =  hdx ,  and 
(13)  simplifies: 
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F  =  2-rro  hu>  /|y!dx  (14) 

iht  °  bc 


The  integral  appearing  in  (14)  is  the  area  between  the 
middle  surface  of  the  annular  section  of  the  shell  and  the 
axis  of  instantaneous  rotation,  as  noted  in  [14] . 

The  energy  dissipated  in  any  of  the  two  plastic  hinges 


is 


h 


h 


•  disc 
Eint 


=  27ia 


/  aoi0dy  =  2iraU0o  /  |y|dy 


(11) 


2  2 

where  a  denotes  the  position  of  the  hinge  with  respect  to  the 

axis  of  symmetry,  and  u@  is  in-plane  component  of  the  velocity  field. 

One  can  easily  evaluate  (15)  as  a  function  of  the  position  of 

the  center  of  instantaneous  rotation  with  respect  to  the  middle 

surface  of  the  hinge  £ : 


£  ^  for  £  =  0  (16a) 

£2  +  |2>  u  <  £  <  |  (16b) 

£h  1  for  £  >  ^  (16c) 

In  the  present  approach  the  hinge  circle  is  indeed  a 

generalized  plastic  hinge  analysed  for  example 

by  Jones  [27]  ,  since  either  bending  or  extension  or  both  can 

be  devleoped  in  it. 


The  total  dissipated  energy  is  the  sum  of  energies 
dissipated  in  continuous  and  discontinuous  dissipation  fields. 
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4 .  Determination  of  the  Crushing  Force 

In  the  actual  solution  the  rate  of  change  of  external  work 
is  always  equal  to  the  internal  dissipation  energy.  A  similar 
equality  for  any  kinematically  admissible  velocity  and  dis¬ 
placement  field  can  be  understood  as  a  definition  of  the  approx¬ 
imate  value  of  the  external  loading.  The  more  close  the 
assumed  fields  are  to  the  real  one,  the  better  the  approximate 
solution  is.  A  considerable  success  of  this  m.ethcd  in  predicting 
the  crush  behavior  of  various  thin-walled  structures  (see  [7,15]) 
is  due  to  the  correct  choice  of  the  admissible  field  with  few 
free  parameters,  and  also  to  a  somehow  weak  dependence  of  the 
solution  on  the  details  of  the  assumed  share  of  the  shell. 

For  example,  any  shape  of  the  shell  middle  surface  between 
inner  and  outer  hinge  circles,  satisfying  kinematic  continuity 
and  giving  the  same  area  in  (14) ,  would  yield  the  same  force 
level.  In  the  present  paper  we  shall  take  the  shape  of  the  shell 
middle  section  in  the  form  of  a  circular  or  parabolic  arc. 
Furthermore,  we  shall  assume  the  deformed  but  rigid  portion  CD 
to  be  a  shifted  mirror  reflection  of  the  original  shape  AB 
(doted  line  on  Fig.  2)  by  a  magnitude  b.  This  implies  opposite 
but  equal  in  magnitude  radius  of  curvature  and  slopes  at  points 
B  and  C.  From  the  continuity  of  displacements  and  slopes  at 
B  and  C,  it  follows  that  the  transition  shape  is  determined  to  within 
a  single  parameter  x>/hich  in  either  the  distance  between  points 
B  and  C,  or  local  radius  of  curvature  r.  Another  free  parameter 
is  the  position  of  the  axis  of  rotation 
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The  rate  of  work  of  external  loading  is 

13  * 

P  -k—  uj  -  for  crushing  between  plates  (17a) 

Pbw  -  for  point  load  or  boss  loading  (17b) 

Now  for  any  given  value  of  the  parameter  describing  the 
crushing  distance  (for  example  central  deflection  or  central 
radius  of  the  annular  zone) ,  the  integral  (14)  is  evaluated  and 

the  value  of  the  force  P  becomes  a  function  of  two  free  para¬ 

meters  £  and  b. 

The  function  P(  b,£)  was  studied  and  it  was  found  that  for  a 
certain  value  of  t,  the  force  attains  a  minimum.  Cn  the  other 
hand,  except  at  a  very  early  stage  of  the  crushing  process,  P 
attains  a  lower  limit  rather  than  a  minimum  with  respect  to 
£  =  0,  Fig.  5. 

Assuming  the  minimum  principle  to  hold  for  geometrically 
nonlinear  problems,  the  lowest  possible  value  of  P  is  identified 
as  the  best  approximate  solution.  This  approach  is  equivalent 
to  the  sequence  of  upper  bound  numerical  calculations  performed 
by  Morris  and  Calladine  [17]  . 

5 .  Tube  Inversion 

The  approximate  solution  of  the  stationary  plastic  process 
of  tube  inversion  with  various  degrees  of  approximation  can  be 
found  in  [11]  and  [15]  .  We  shall  study  the  same  problem  to  explain 
the  details  of  the  derivation  and  to  determine  the  accuracy  of 
the  present  approach  by  comparing  it  with  known  results. 
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Consider  a  tube  of  the  initial  radius  R  being  inverted 
inside  out  to  a  final  radius  R  +  2r.  At  the  hinge  circle  B 
the  straight  line  generator  of  the  cylinder  is  bent  into  the 
curvature  p  and  at  the  hinge  circle  C  the  curvature  is  re¬ 
moved,  Fig.  6.  A  simple  and  yet  accurate  solution  of  this 
problem,  accounting  for  continuous  circumferential  extension  and 
concentrated  bending  at  B  and  C,  was  worked  out  by  Kitching  et  al .  [15] 
and  has  the  form 


Pmin  .  1 3r  ~  ,  QO  fR.  fP-h 

2rtT  *  4  Vt  ’  6-92  Vh'  r  "VT  <1S) 

In  order  to  apply  the  present  approach,  we  assume  that  the 
portion  of  the  shell  between  hinge  circles  is  and  C  form  one 
naif  of  the  toroidal  surface.  We  also  assure  for  simplicity 
chat  r/h  >  1,  so  that  Eq.  (14)  applies  and  t,  -•  0.  The  distance 
between  B  and  c  is  h  =  2r  and  the  integral  appearing  in  (16) 
represents  the  area  of  the  half  circle.  Using  this  result  and 
Eq.  (15),  the  rate  of  energy  balance  equation  yields 

•2  ’  .h2 

Pr  oi  =  2t to  hw—s—  +  2tto_uA-  [R  +  (R+2r)  ]  (K  ) 

O  2.  0  4 

The  rate  of  infinitesimal  rotation  w  drops  out  from  both  sides 
of  (19)  and  the  dimensionless  force  becomes 


(?-C) 


f’inimizing  the  right  hand  side  of  (20)  with  respect  tc  r  and 
substituting  the  optimal  radius  back  into  (2C)  we  can  obtain 


the  final  formula: 
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j2i2  -  4^2  +  2  =  7.09^  ♦  2.-  r  =Vt  ,21> 

The  above  results  compares  well  with  solution  (18) ,  especially 
for  large  R/h. 

It  is  interesting  to  estimate  the  effect  of  the  variable 
circumferential  curvature,  according  to  the  approximate  formula 
(11) .  One  can  easily  find  that  the  magnitude  of  the  non- 
dimensional  force  in  (21)  increases  then  by  2,  which  again  is  a 
small  contribution  for  large  R/h.  At  the  same  time  the  optimum 
radius  of  inversion  remains  the  same. 

It  can  be  checked  that  the  condition  (10)  of  kinematic 

V  V 
B  C 

continuity,  reduces  now  to  —  =  —  and  thus  is  satisfied. 

Finally  it  should  be  noted  that  the  crushing  force  of  a 
tube  loaded  by  a  central  force  rather  than  crushed  between 
plates,  is  one  half  of  that  predicted  by  (21). 

6 .  Crushing  of  a  conical  shell 

The  geometry  of  the  partially  crushed  conical  shell  is  shc'x 

in  Fig.  7.  We  assume  that  the  portion  of  the  shell  subjected  to 

localized  plastic  deformation  forms  a  section  of  a  toroidal  curface 

and  that  the  inverted  cone  is  of  the  same  angle  a.  The  area  of 

r2 

the  section  of  a  circle  between  the  points  E;  and  C  is  S  =  y  [(jr-2a) 
-  sin(7r-2a)]  and  the  distance  b  =  2rcosa.  It  is  now  straight¬ 
forward  to  work  out  the  solution  which  takes  the  form 

P  — 

Vh'  V(7T"2a)  '  sin(7T-20l)  (22) 
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-V 


ha 


(iT-2a)  -  sin(ir-2a)  (23) 

Here  a  is  the  measure  of  the  crushing  distance,  which  is  related 
to  the  central  deflection  by 


w  = 


2a 

tana 


(24) 


The  formulas  become  increasingly  inaccurate  for  a  -r-  and 
also  for  small  since  then  the  assumptions  leading  to  5  =  0 
are  not  satisfied.  In  these  cases  a  more  exact  solution  can 
be  obtained  when  needed,  using  (13)  rather  than  (14)  and 
minimizing  the  solution  with  respect  to 

Equation  (23)  predicts  an  increasing  radius  of  the  toroidal 
surface,  r  >  0,  as  the  crushing  zone  is  expanding.  The  change  of 
the  radius  is  related  to  velocities  of  hinge  circles  by 


2r  =  (VB~Vc)tana 


(2i) 


At  the  same  time  we  have  assumed  that  =  Rc  =  »  and  rg  =  r^ 
which  yields  V  =  V  .  We  see  that  the  condition  of  kinematic 

D  L. 

continuity  is  not  satisfied  by  the  present  solution.  The 
difficulty  encountered  can  be  overcome  by  considering  the  general 
case  in  which  r_  f  r_  R_  ^  ®.  Now,  the  toroidal  surface  ceases 

D  v-  C 

to  be  of  a  circular  cross-section  and  the  hinge  circles  B  and  C 
will  no  longer  be  positioned  on  the  same  level.  As  a  result 
membrane  force  will  develop  in  the  hinge  C  which  would  result  in 
a  slight  increase  of  the  force  level.  Also,  with  a  variable 
radius  R  (a)  the  hinge  circle  will  produce  a  curved  line  rather  than  a 
conical  one  on  the  meridional  olane.  It  is  not  difficult  to  write 
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a  system  of  ordinary  differential  equations  with  a  shifted 
argument  describing  variation  of  Vn  V_  r_  r_  and  R_  with 
the  crush  distance.  We  shall  not  attempt  to  study  this  sytem 
since  it  is  believed  that  insiginif icant  changes  would  have  then 
been  introduced  to  the  solution.  The  problem  however  is 
interesting  by  its  own  right. 
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Equations  (26)  and  (29)  furnish  a  parametric  representation  of 

the  function  force  versus  displacement.  A  plot  of  this  function 

for  several  values  of  the  radius  to  thickness  ration  R/h  is  shown 

in  Fig.  9.  The  dependence  of  the  solution  on  R/h  appears  to  be 

h 

weak.  Indeed,  treating  ^  as  a  small  parameter  and  expanding  the  square 

h  h 

root  in  (29)  in  power  series  with  respect  to  —  arcund  —  =  0  one  gets 

2  2 

1  -  4!  — - 1  -  2  R  +  •  •  •  (30) 

so  that  (29)  can  now  be  approximated  by 

_  .  2 

w  ^  sin  ct  /n1N 

h  2K  a-cosa)  -  2  — (31) 

Combining  next  (26)  and  (31)  to  eliminate  the  parameter  R/h 
we  obtain 

w  _  /pmin  \  1 1  -ccso-l  n  (  sin?'o) 
h  \2Mott+7  |2n(a)j  z  |  n  ( a )  j 

Expanding  the  trigonometric  function  appearing  in  (32)  in  power1 
series  and  retaining  terms  up  to  the  third  power,  the  two 
coefficients  in  (32)  become 


1-cosg  _  _1  2sina  -  sin2g 
2n  (a)  4  2a  -  sin2a 


2sin3a  2s 


in3a  a  3  (,  a A 
-  sin2a  2  \  6  / 


(33a) 


(33b) 


Thus,  the  first  coefficient  turns  to  be  constant,  while  the 
second  varies  very  little  which  a,  especially  for  a  <  1.  Sub¬ 


stituting  (33)  back  into  (32)  ,  the  final  formula  for  the  force- 
deflection  relationship  takes  the  form 


The  plot  of  the  function  (34)  is  given  in  Fig.  9  by  a  broken 
line.  The  agreement  with  the  exact  solution  (full  line),  which 
depends  on  R/h,  is  seen  to  be  very  good. 

It  is  interesting  enough  that  the  same  solution  can  be 
obtained  by  making  a  different  set  of  assumptions.  Suppose  the 
hemispherical  shell  is  approximated  by  a  parabolical  shell.  The 
relationship  between  the  angle  a  and  central  deflection  is  now 
given  by 


w 

h 


R 

h 


2 


sin  a 


(35) 


Replacing  likewise  the  circular  cross-section  of  the  toroidal 
surface  by  a  parabolic  arc 


the  area  integral  in  (14)  can  be  readily  evaluated,  and  the 
expression  for  the  crushing  force  takes  the  form: 


(36) 


2^rM, 


2b  2R  sina  i 

-  3  h  Slna  +  - E - 1 


(37) 


Optimization  the  right  hand  side  of  (37)  with  respect  to  b 
’-esults  in 


b  =  3hR 


(38) 


Now  combining  (35),  (37)  and  (38),  the  expression  (34)  for  the 

crushing  force  is  obtained  in  the  same  form  as  before.  This  result  is 
important  from  the  practical  point  of  view,  since  the  mathematics 
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involved  in  the  alternative  derivation  of  (34)  are  much  simpler. 
We  shall  take  advantage  of  this  property,  and  in  the  remainder  of 
the  paper  dealing  with  spherical  shells  always  replace  a  circle 
by  a  parabolic  arc.  Note  that  by  doing  so  the  condition  of  con¬ 
tinuity  of  slope  of  the  displacement  field  (6)  is  not  satisfied, 
but  as  we  have  seen,  the  solution  is  hardly  affected  by  this. 

The  second  continuity  condition  (10)  reduces  now  to 


V 

V, 


B 


ill 

*♦5 


(39) 


which  means  that  |  V  I  f-  1 V  j  .  On  the  other  hand  formula  (38) 

b  L 

in  the  present  solution  implied  that  b  f  0.  Recalling  that 
b  =  ( V  -  V_ )  cosa,  we  have  I V  |  =  | V„ ' ,  which  contradicts  the 
condition- of  continuity.  A  consistent  set  of  equations  can 
be  written,  as  discussed  previously.  However,  the  resulting 
solution,  which  apparently  looses  its  appealing  simplicity, 
is  believed  to  introduce  insignificant  changes  and  thus  will  not 
be  discussed  further. 

The  present  solution  with  two  moving  hinges  becomes  increas¬ 
ingly  inaccurate  with  decreasing  value  of  the  parameter  Two 

factors  are  responsible  for  that.  First,  at  an  early  stage  of 
the  crushing  process,  only  one  hinge  circle  is  formed  and  the  kine¬ 
matics  of  the  process  does  not  involve  a  free  parameter  to  be  opti¬ 
mized.  Secondly,  for  small  ot  the  portion  of  the  toroidal  shell 
becomes  relatively  shallow  so  that  the  assumption  leading  to 
l  -  0  is  no  longer  valid,  and  the  shell  should  be  treated  as  a 
three-dimensional  body,  according  to  the  method  developed  by 
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Calladine  [18].  Te  shall  study  these  two  cases  separately. 

Assume  that  the  radius  of  the  inner  hinge  is  zero.  From 
simple  geometrical  consideration  Fig.  10/  it  follows  that 

b  =  R  sina  ,  R  =  2  r 

Introducing  (40)  into  (35)  and  (37)  and  eliminating  R/h  we 
obtain  the  following  force-deflection  characteristics 


4_  w 
3  h 


1 


which  is  identical  to  the  one  derived  in  [27] •  The  function 


(40) 


(41) 


(41)  represents  a  straight  line  intersecting  the  curve  described 
Ly  (34)  at  g  =  |,  Fig.  11. 

Return  now  to  the  general  case  of  the  shell  in  which  the 


internal  energy  dissipation  is  governed  by  (13)  and  (16b) .  We 
assume  that  (40)  holds,  and  that  the  shape  of  the  shell  middle 
surface  represents  a  parabolic  arc.  The  energy  balance  equation 
can  be  written  in  the  form 


b 


■•here  yl  and  yll  are  respectively  the  upper  and  lower  contour  of 
the  shell  cross-section,  measured  from  the  axis  of  instantaneous 
rotation,  see  Fig.  4. 


R 


(43a) 
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y 


ii 


(43b) 


The  integrals  in  (42)  can  be  readily  evaluated  and  the 
crushing  force  becomes  a  function  of  a  single  parameter  £ 


P 

2ttM 

o 


8 


15  R  h 


Introduction  (40)  into  (35)  we  obtain  a  unique  relation 
between  the  radius  of  the  hinge  and  central  deflection 


,2 

w  _  b 
h  2hR 

The  minimization  of  the  right  hand  side  of  (44)  yields  an 
optimum  value  of  the  parameter  £ 


£  _1_  b^  1  w 
h  ~  12  Rh  6  h 


(44) 


(45) 


(46) 


Substituting  (45)  and  (46)  back  into  (44)  we  finally  obtain 


P 


2ttM 

o 


2 


(47) 


The  above  relation  is  depicted  in  Fig.  11  by  a  solid  line. 


-25- 


Ecuation  (47)  is  valid  as  long  as  the  line  of  instantaneous 
rotation  first  touches  either  the  inner  or  cuter  contour 
y^  cr  v^  .  This  first  occurs  at  ^  Thus,  the 

range  of  applicability  of  both  solutions  appears  to  be  the  same. 
The  formula  (41)  is  easy  to  evaluate  but  it  gives  unsatisf actoy 
results  for  w/h  =  0.  This  will  become  even  more  evident  in  the 
problem  of  a  shell  loaded  ,’y  a  rigid  boss.  Cn  the  other  hand,  the 
solution  (47)  predicts  correctly  the  yield  point  load  but  the 
mathematics  involved  are  much  more  complicated. 

The  4%  difference  between  the  functions  (34)  and  (47)  at 
w/h  =  —  is  attributed  to  the  fact  that  the  former  involves  the 
simplifying  assumption  z,  =  0 ,  while  the  latter  not.  By  con¬ 
tinuing  the  exact  calculation  beyond  that  point,  the  corresponding 
solution  (doted  line  in  Fig.  11)  would  gradually  converge  to  the 
curve  representing  the  formula  (34 ) .  The  thin  full  line  in  the 
same  figure  denotes  the  Morris  and  Calladine  numerical  solution. 
Unfortunately,  no  information  was  given  in  [17]  on  the  magnitude 
of  the  integration  step  employed  so  it  is  difficult  to  assess 
its  accuracy.  A  good  agreement  of  both  solutions  is  observed 
over  the  entire  range  of  deflection  covered  by  the  Morris  and 
Calladine  solution. 

t .  Crushing  of  a  Spherical  Shell  Between  Rigid  Plates 

The  procedure  for  solving  this  problem  essentially  follows 
that  presented  in  the  preceding  section,  except  that  now  the 
rate  of  external  work  should  be  computed  according  to  (17b)  and 
the  relation  between  w  and  a  is  different.  As  before, 
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the  analysis  follows  several  stages.  The  yield  point  and  post¬ 
yield  behavior  is  accurately  predicted  by  the  formula 


=  4  + 


!  (h) 


w  3 

valid  for  ^  ^  .  In  the  later  stage  of  the  crushing  process, 

when  the  inner  hinge  leaves  the  shell  axis,  the  plate  displacement 
is  related  to  the  angle  a  by 


w  R  .  2  h 
_  =  _  sin  a  _  _  sina 

The  optimum  value  of  width  of  the  deforming  radius  is  the  same 
as  before  b  =  /3hR,  Equation  (38),  and  the  resulting  force- 
deflection  relationship  takes  the  form 

w  _  _1_  /  P  .  ,\2  _3_  /  P  . 

h  “  43  \2ttM0  ■  )  "32  \2ttMq 

The  above  solution  is  valid  theoretically  until  the  crush 

distance  reaches  the  shell  radius  w  =  R.  It  is  interesting  to 

compare  the  formula  (50)  with  an  approximate  solution  due  to 

Updike  [14]: 


which  was  believed  by  the  author  to  be  valid  for  deflections 
not  exceeding  about  one  tenth  of  the  shell  radius.  When  both 
solutions  are  plotted,  it  is  seen  that  a  relative  difference 
between  them  is  small  over  the  entire  range  of  variation  of  w, 
Fig.  12. 


(49) 


(50) 


(51) 
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9 .  Rigid  Boss  Loading 

Leckie  and  Penny  performed  a  series  of  tests  on  carefully 
manufactured  spherical  shells  with  R/’n  =  60  and  variable  boss  size 
[16].  They  observed  that  with  increasing  downward  deflection 
the  crushing  force  for  a  given  radius  of  the  boss  p  reaches  a 
maximum,  then  drops  suddenly  as  the  shell  buckles,  and  starts 
to  rise  again.  They  also  found  that  the  initial  plastic  yield- 
point  load  depends  markedly  on  a  single  geometrical  parameter 

B  =  (52 

/Rh 

but  with  increasing  deflection  this  dependence  is  much  less 
pronounced.  However,  the  tests  were  not  run  far  enough  to  decide 
whether  all  P  -  w  curves  did  converge  to  a  single  one. 

Assuming  a  rigid  perfectly  plastic  material  idealization 

Morris  and  Calladine  [17]  were  able  to  predict  accurately  the 

initial  response  of  the  shell.  Their  numerical  calculations  have 

shown  that  while  the  slope  of  the  load-deflection  curves  becomes 
w 

the  same  for  ^  _>  2,  the  curves  corresponding  to  different 
3  are  considerably  shifted  with  respect  to  one  another. 

It  is  reasonable  to  expect  that  the  effect  of  the  boss 
size  is  a  local  one  and  should  be  almost  entirely  "forgotten" 
by  the  shell  at  later  stages  of  the  deformation  process,  so  that 
all  curves  would  eventually  converge.  Indeed,  the  deformation 
mechanism  with  expanding  ring  of  plastic  deformation  is  the  same 
as  before,  except  that  the  reference  point  for  measuring  deflec¬ 
tion  changes  with  the  boss  size.  This  gives  rise  to  a  new  term 
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in  the  force-displacement  relationship 

2 

+  i: 


X  =  J_  /_L_  +  A  -  3  . 
h  16  \2ttM  +  2  2 


(53) 


The  above  solution  reduces  to  the  one  describing  the  shell 
under  central-point  load  (formula  (34))  by  setting  £  =  0.  The 
new  term  is  relatively  small,  and  provides  a  shift  of  P  -  w 
curves  in  the  w/h  direction  Fig.  13.  The  horizontal  displacement  of 
P-v;  curves  in  the  numerical  solution  reported  in  [17]  for  the  same 
values  of  8  are  much  larger,  which  is  probably  attributed  to 
the  trial  and  error  optimization  procedure  employed . 

We  were  unable  to  derive  a  simple  closed-form  solution  for 
the  initial  phase  of  loading  based  on  the  exact  formula  (13) 
showing  the  nature  of  the  convergence.  However,  the  interesting 
pattern  of  behavidr  of  the  shell  associated  with  plastic  snap- 
through  can  be  described,  at  least  qualitatively  using  the 
approximate  formula  (14) .  Assuming  that  a  stationary  hinge  circle 
is  formed  at  T  =  p  and  the  moving  one  at  T  =  p  +  b,  the  load- 
deflection  relationship  was  found  to  be 


P 

2itM 


-B 


Ft* 


2  B 


FT 


+  6-8 


1 


(54) 


The  above  solution  reduces  to  (41)  for  8=0  but  it  leads  to  an 
unrealistic  infinite  force  magnitude  with  ^  0 .  This  is  a 

consequence  of  the  assumed  symmetric  geometrical  shape  of  the 
shell  cross-section  (36)  with  zero  initial  width  of  the  plastic 
zone  to  initiate  the  motion.  In  reality,  the  shape  is  unsymmetric 
with  finite  width.  It  was  possible  to  derive  an  approximate 
expression  for  the  initial  yield-point  load  (w/h  =  0)  in  the  form: 
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2ttM 


=  86,  b  = 


/Rh 


Rh 

2p 


(55) 


The  plot  of  the  function  (54)  with  (55)  for  several  values  of 
the  parameter  6  is  shown  in  Fig.  14. 

10 .  Spherical  Cap  Under  External  Pressure  Loading 

Consider  a  simply  supported  spherical  shell  with  base  radius 
a,  and  rise  H,  Fig.  15.  The  shallowness  of  the  shell  is  described 
by  the  parameter 


a 


(56) 


Following  Ref.  [21]  define  the  total  load  P  and  the  reference 
load  PQ  by 

P  =  -rpa2  P  =  6ttM  (57) 

c  o  o 

where  p  is  a  uniformly  distributed  downward  pressure.  It  is 
relatively  easy  to  derive  the  load-deflection  relationship  for 
any  value  of  a,  using  the  present  method.  We  shall  however 
restrict  the  analysis  to  the  limiting  cases  of  small  and  large 
a.  Consider  first  the  case  a  >>  1  for  which  equation  (14)  applies. 
Suppose  the  outer  boundary  of  the  plastically  deforming  zone  BC 
is  b  while  the  portion  AB  remains  rigid.  The  assumed  velocity 
field  consist,  as  previously,  of  the  rotation  around  the  point  B 


•  • 

w(x)  =  tux  (x  1  b  £  a)  (58) 

This  rotation  imposes  a  continuous  change  in  the  initially 

circular  shape  of  the  shell  middle  surface.  If  the  central 

deflection  is  smaller  than  the  rise  of  the  shell  above  the 
b2 

circle  b,  w  <  ,  the  current  shape  can  be  approximated  by 
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the  parabolic  arc. 

y(x)  =  (l)~  (Ir  +  w) 

For  the  time  being  we  assume  the  plastic  hinge  circle  to  be 
located  somewhere  between  the  shell  center  and  the  support 
0  <  b  <_  a.  Equating  the  rate  of  external  work  to  the  energy- 
dissipated  in  the  region  of  continuous  deformation  and  plastic 
hinge  circle  one  gets 


2  7T 


/ 


pw(x)xdx  =  2irq 


70ha)  /  y(x 


)  dx  +  2 ttoq  —  mb 


or  after  integration 


16  0) 


P_ 

P 

o 


i  Si  +  i 

3  h  +  1 


(61) 


Id 

where  b  =  JhR'  We  can  see  for  t*ie  force  to  be  a  minimum 

one  has  to  consider  the  largest  possible  value  of  b,  i.e.  b  =  a. 

With  the  assumed  simply  supported  boundary  conditions  there  will 
not  be  any  contribution  of  the  concentrated  hinge  to  the  total 
dissipation  energy,  and  the  last  term  in  Eq.  (61)  should  vanish, 
giving  the  final  expression  for  the  force  deflection  characteristics 

—  =  ja  -  for  w  <  a  (6  2) 

o 

In  particular,  the  initial  yield-point  load  is  obtained  from 
(62)  by  putting  w  =  0,  P/P0  =  8/3a.  A  plot  of  the  above  formula 
is  undistinguishable  from  the  curve  representing  the  much  more 
complicated  solution  derived  by  Kondo  and  Pian  [21],  Fig.  16 

At  the  other  extreme  we  have  very  shallow  shells  in  which 
the  line  of  instantaneous  rotation  lies  within  the  contour  of 
the  shell  cross-section,  Fig.  4.  This  implies  a  <  1. 
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Assume  that  the  shape  of  the  shell  middle  surface  is 
changing  with  deflections  according  to  Eq.  (59).  The  energy 
in  the  continuous  deformation  field  is  given  by 


E  ,  =  2ira_ 

cont.  o 


fh\  ryI  / 

~yZI  ) 

/  J  J 

ydy  j 

where  the  limits  of  the  integrals  appearing  in  (63)  are 

I  h  r 

y  =  Y  +  2  “  ? 


II  ^  h  ^  r 

y  =  -y  +  ,  +  C 


The  expression  for  the  energy  dissipated  at  the  outer  hinge 
would  depend  on  the  position  of  the  hinge  circle.  If  b  <  a, 
then  the  boundary  conditions  do  not  have  any  effect  on  the 
solution  and  Eq .  (16b)  applies.  If  the  plastic  zone  extends 

up  to  the  support,  b  =  a,  then  for  a  simply  supported  shell 
restrained  from  axial  motion,  the  term  responsible  for  bending 
vanishes 


E,  .  =  2TTG'_tok 

hinge  c 


rh2  + 


n  = 


J 


1  if  b  <  a 
( 0  if  b  =  a 

The  left  hand  side  of  the  energy  balance  equation  is  the  same 
as  before  (60) ,  while  the  right  hand  side  is  a  sum  of  (63)  and 
(65).  Evaluating  the  integral,  one  obtains 


P  b 

P  a 
o 


(1  +  H)  +  ||b2-^bq  +  |-q2  +  2C2-:j5(2b-q) 


where 


*  q  = 


w 

h 


(63) 

(64a) 

(64b) 


(65) 


(66) 
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After  minimization  with  respect  to  £,  the  formula  (66)  becomes 


P  b  ,  56  r2  68  t—  26  2 

p  a  ~  (1  +  n)  +  45  b  45  bq  45  q 

o 

The  load-deflection  curve,  represented  by  (67)  is  dropping 
with  increasing  b  and  attains  the  lowest  position  when 


(67) 


-2  45  A  26  , 

b  =  (1  +  n)  56  +  5^  g2 


(68) 


For  example  at  q  =  0,  the  optimum  value  of  radius  of  the  outer 
hinge  is  b  =  1.3,  falling  beyond  the  range  of  applicability  of 
the  present  solution,  which  is  b  <  a  <  1.  Consequently,  the  least 
value  of  the  crushing  force  is  obtained  by  taking  the  largest 
possible  radius,  i.e.  b  =  a.  We  have  thus  shown  that  for  shallow 
shells  the  zone  of  continuous  plastic  deformations  is  extending 
right  up  to  the  support.  Substituting  n  =  0  and  a  =  b  into  (59) , 
leads  to  the  following  force-deflection  relation. 


_P 

P 


56  2 
45 


—  =  1  +  —  OL  -  Tr-  Ctq  +  —  q 


68 

45 


26  _2 
45 


(69) 


In  particular,  the  initial  yield-point  load  is  given  by 

P  ,  .  56  2  . 

P“=i45a  a<l 

o 


(70) 


In  deriving  Eq.  (69) ,  shell  deformations  both  in  radial  and 
circumferential  directions  were  taken  into  account.  A  relative 
contribution  of  the  two  components  of  the  dissipated  energy  is 
controlled  by  the  position  of  the  axis  of  instantaneous  rotation  I~I . 
The  optimum  position  is  one  which  minimizes  the  magnitude  of 
the  pressure  required  to  maintain  the  plastic  flow.  It  is 
interesting  to  note  that  the  neglect  of  meridional  strain,  an 
assumption  made  in  Ref.  [21],  would  lead  tc  an  increase  rather 
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than  decrease  in  the  level  of  the  yield-point  load.  Indeed,  in 
the  present  approach  all  work  done  by  the  meridional  strain  is 
concentrated  in  the  generalized  tensile  hinges  and  the  corres¬ 
ponding  dissipation  energy  is  given  by  (65) .  This  term  can 
easily  be  made  equal  to  zero  by  taking  5=0.  Substituting 
a  =  b  and  5=0  into  (66)  we  obtain  a  new  solution 


P 

P 


.  x  32  r2  12  ,4  2 

=  1  +  —  b  _  -F-  aq+^-q 


(71) 


o  -  15  “5 

Load-deflection  curves  for  few  chosen  values  of  the  parameter  a 
are  shown  in  Fig. 17  together  with  the  solution  due  to  Kondo  and  Pian. 
The  initial  yield-point  load  is  obtained  from  (71)  by  setting 


q  =  0  : 


p  i  32  2  ^ 

=  1  +  rj-jr  a  a  £  1  <72) 

o 

The  above  formula  is  identical  to  the  one  derived  by  Kondo 
and  Pian  and  the  corresponding  P  -  a  curve  clearly  lies  above 
that  given  by  (70),  Fig.  16.  Thus,  the  neglect  of  meridional 
strain  rate  was  offset  by  even  larger  increase  of  the  dissipa¬ 
tion  due  to  the  circumferential  strain  rate.  With  increasing 
ct,  the  contribution  of  dissipated  energy  in  meridional  direction 
diminishes,  illustrating  a  nature  of  the  approximation  introduced 
in  replacing  (13)  by  (14).  Note  that  the  solution  (62)  was  derived 
on  the  assumption  that  5  =  and  this  explains  the  excellent 
correlation  of  the  present  solution  with  the  more  elaborate 
solution  due  to  Kondo  and  Pian,  which  employs  a  full  set  of  equations 
describing  moderately  large  deflection  of  shells. 
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It  is  not  in  principle  difficult  to  extend  the  present 
solution  to  shells  with  intermediate  values  of  the  parameter  a. 
The  P  -  a  curve  would  first  follow  Eq.  (70)  and  then  converge 
assymptotically  to  Eq .  (62). 


11 .  Effect  of  End  Condition 

According  to  the  present  analysis,  the  boundary  conditions 
are  not  felt  by  the  shell  until  the  outer  hinge  reaches  the 
support  located  at  p*.  It  is  seen  from  (35)  and  (38)  that  this 
occurs  when  deflections  become 


(73) 


For  a  shell  under • central-point  load,  the  width  of  the  plastic 

zone  predicted  by  the  approximate  solution  is  b  =  /3Rh .  From  now 

on  the  outer  hinge  remains  fixed,  and  any  further  increase  of  w/h 

is  due  to  the  diminishing  of  the  width  of  the  toroidal  section  b. 

Introducing  the  dimensionless  parameter  a  =  — ^  ■ ,  the  geometrical 

/3Rh 

relation  becomes 


w_p*2  32  w*  3 

h  ~  TT  ~  2  a  =  TT  +  7 


(1 


2 


a 


) 


For  fully  clamped  boundary  conditions,  the  solution  (37)  with 
two  concentrated  plastic  hinges  still  applies,  so  that  after 
introducing  (73)  into  (37)  we  have 


P 

2ttM0 


/3Rh 


+  b 


i 


(74) 


(75) 


Equations  (74)  and  (75)  provide  a  parametric  representation  of  the 
load-deflection  relationship  at  the  terminal  stage  of  the  shell  motion 
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Kith  a  -*■  C  the  force  blows  up  to  infinity  . Fig-  18.  This  occurs 
when  +  j.  In  reality  the  force  level  in  shells 

subjected  to  central  point  load  is  equal  to  the  strength  of  the 
inverted  spherical  shell  in  tension.  An  upper  bound  for  this 
force  is 


Ptens  _  4R 
2  Mqt  h 


(76) 


This  is  much  higher  than  (75)  with  a  -  1,  but  finite.  For  shells 
compressed  between  rigid  plates  the  formula  (75)  applies 
multiplied  by  the  coefficient  2.  The  force  increases  with  a 
approaching  zero  and  indeed  goes  to  infinity  when  the  upper  plate 
comes  in  contact  with  the  supporting  plate.  A  typical  loading 
situation  for  the  two  cases  described  above  is  shown  in  Fig.  18. 

12 .  Conclusions 

This  paper  presents  a  general  methodology  for  predicting  the 
crush  resistance  of  arbitrary  rigid  plastic  shells  of  revolution 
deforming  in  an  axisymmetric  mode.  The  implications  of  the  condi¬ 
tions  of  kinematic  continuity,  as  well  as  the  boundary  conditions, 
are  discussed.  Several  practical  problems  are  solved,  including 
a  conical  shell  and  a  spherical  shell  crushed  between  rigid  plates 
and  under  boss  loading,  and  a  spherical  cap  under  external  uniform 
pressure . 

For  the  initial  stage  of  deformation  the  structure  has  been 
treated  as  a  three  dimensional  body,  in  the  spirit  of  Calladine's 
approach  [18],  while  in  subsequent  stages  a  simplification  to  this 
approach  has  been  suggested,  which  results  in  simple  closed  form 
solutions  useful  for  engineering  applications.  A  good  correlation 
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was  found  with  some  existing  solutions  of  similar  problems  obtained 
through  much  lengthier  calculations. 

In  all  cases  treated,  the  zone  of  plastic  deformation  was 
found  to  be  confined  to  a  relatively  narrow  region.  Also,  the 
level  of  the  crushing  force  was  shown  to  be  quite  insensitive  to 
the  particular  shape  of  that  2one . 

The  present  method  makes  it  possible  to  study  the  crush 
behavior  of  rotationally  symmetric  shells  strengthened  by  circum¬ 
ferential  ring  stiffeners.  The  stiffner  acts  initially  as  a 
clamping  ring  until  the  force  reaches  a  critical  level.  Then  the 
ring  becomes  a  rigid  boss.  Thus  the  determination  of  the  initial- 
point  load  for  a  shell  loaded  by  a  rigid  boss  provides  a  basis 
for  studying  the  optimum  design  of  stiffened  plastic  shells 
*  against  crash.  In  the  optimum  design  the  rigidity  of  the  stiffeners 
must  be  compatible  with  the  local  crushing  strength  of  the  shell. 
This  problem  will  be  studied  in  detail  in  a  future  publication. 
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